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1. Introduction 


As one of the theoretical possibilities for the initiation of frontal 
rainbands, SI (Symmetric Instability), moist SI, and especially CSI 
(Conditional SI) have been shown to have many attractive features (Hoskins, 
1974; Raymond, 1978; Bennetts and Hoskins, 1979; Emanuel, 1979, 1983; Parsons 
and Hobbs, 1983, Xu and Zhou, 1982; Xu, 1984). However in comparison with BI 
(Buoyancy or Convective Instability) (cf. Chandraseckhar, 1961 and Joseph, 
1976a, b), the SI problem has not been explored sufficiently. Obviously the 
presence of basic wind shear and, by the thermal wind relationship, a 
horizontal temperature gradient makes SI problems generally more difficult 
than BI problems. 

It has been noted, however, that the most unstable motion of SI may be 
locally described as II (Inertial Instability) on isentropic surfaces in an 
hydrostatic atmosphere (Hoskins, 1974) and that a linear viscous SI problem 
can be transformed into a form of linear viscous BI or II (Xu and Zhou, 1982). 
Thus the similarity and intrinsic relationship among SI, BI and II for 
nonlinear viscous flows becomes more perceptible. In the following section we 
will see that the governing SI equations can be changed via a coordinate 
transformation into a form similar to those of the BI problem. Thus the 
well-studied BI theory can be utilized to improve our understanding of SI 
problems. 

Since a similarity among SI, BI and II exists, it will be expected that 
SI has the same basic nature of BI and II and thus the physical picture of SI 
may be illustrated almost as intuitively as that of BI or II. In fact, by 
treating buoyancy and inertial effects separately, an intuitive analysis for 


the parcel dynamics of SI has been attempted by Emanuel (1983). A similar 
idea was also used by McIntyre (1970) to explain the diffusive des tablizat ion 
mechanism in a viscous SI problem where the Prandtl number differed from 
unity. Referring to Emanuel's analysis, we will emphasize the analogy between 
SI, BI and II. To illustrate the properties of SI, we shall analyze the 
dynamics and energetics of the parcel motions in Section 3. We shall also 
examine to what extent the results for dry SI may be applied to moist SI. 

Moist SI will be discussed in Section 4, where a qualitative analysis for CSI 
will also, be involved. 

2. Similarity Between SI, BI, and II 

Assume that V = (0,V(x,z),0) and 0 Q (x,z) are the basic wind and potential 
temperature fields respectively (cf., Fig. 1). By the Boussinesq approximation 
we have the following nonlinear viscous equations for symmetric perturbations 

7 1 - fv + a ^ = a 7 • p v7 u , (2.1a) 

dt o dx o o 
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where f is the Coriolis parameter (assumed constant) and P Q (x,z) = a ^ (x,z) 

d _ 3 3 3 

is the density of the background atmosphere. Here — = + u + w , 

0 • 

(u,w,v,9,p) = (u' ,w' ,v' ,g — ,p' ) and the primes represent perturbation 

o 

quantities. In (2.1b) the effect of perturbation pressure on buoyancy 
is neglected. The symbol v is the kinematic molecular " or eddy diffusivity. 

N, S and F the Brunt Vaisala, baroclinic and inertial frequencies respectively, 
are defined as follows: 
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where thermal wind relation has been used. 

Consider now the local behaviour of (2.1) by assuming that 
and are constants. We can introduce a rotational transformation 

a * 

_ cos ot in& ^ 

= ( . * A ) where positive a corresponds to anticlockwise rotation, 

ct sinoc cosot 

i. e. , see Fig. 1, 
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where 


A 1 . , -2S N 

a = Y Arctan (— -z j) 

N - F 


(2.4) 


With (2.3) and (2.4), (2.1) can be transformed into 


(— — + v * V - a V • P v7)v - h + a = 0 , 

o t O O ~ ~ O 
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where 7 = (3/3x, 3/3z), y = (u,w), = (a,b), and n = ( is the 

2 .2 0 N 

F S 

diagonal form of stability tensor H = ( ). Here 

S Z N 
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It has been shown by Xu and Zhou (1982) that (2.4), or equivalently 
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gives the most unstable or least stable slope angle for linear SI, i.e., 
a disturbance which induces parcel displacements at this angle in the 
(x,z) plane will be most unstable. Now with (2,5) it is easy to see that 

a 

the most unstable angle a is indeed only a property of the basic flow and 
not changed by including nonlinearity. From (2.5) we can see that if 

< F^, then z is in the direction of a with ■ min u>2(a) and 

-2 2 a 
F = m^x to (a) where 

2 2 2 2 2 2 
to (a) = N sin a + 2S sinot cosa + F cos a (2.8) 

(cf. Ooyama, 1966, Hoskins, 1974, or Xu et al., 1982). If > F^ , then 

~ ^ 2 *2 2 
x is in the direction of ol where N * max w (ot) and F = min w (a). 

a a 

2 2 

In fact, the above two cases, corresponding to the sign of N - F , share a 
similar nature symmetrically because one can change into the other simply by 

AAA AAA 22 2 

interchanging (x,u,a) and (z,w,b). For N < F and F >0, Eqs. (2.5) have 

a similar form to the governing equations of nonlinear BI problems in a 
~2 a 

stratified (N ), rotating (F) fluid. Thus, even in the nonlinear sense, SI 

A <A 2 

is similar to BI and z is the most unstable (if N < 0) or least 

- 2 a 2 2 

stable (if N >0) direction, i.e., the direction of oi. For N > F and 
2 

N >0, Eqs. (2.5) indicate that SI is similar to II (Inertial Instability). 

However, as we mentioned above, these two cases corresponding to the sign 
2 2 

of N - F are similar. Thus, BI, II and SI are similar, at least, if 

A A A A 

we only consider the transverse circulation (u,w) and take (a,b) as two 
scalars. 
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So far we can see that, although BI and II are particular cases of SI, 
the local SI or SI problem with a uniform basic state can be transformed into 
the form of BI or II. Furthermore both linear (Xu and Zhou, 1982) and 
nonlinear theories (cf. Chapter IV, Xu, 1984) of generalized energetics 
indicate that the nonuniformity of the basic state does not change the basic 
nature of BI, II and SI. Therefore, it is indeed feasible either to induce 
their common natures from already known II, SI, and especially BI theory, or 
to deduce certain properties for a particular stability problem, especially 
the SI problem which we are particularly concerned with from their common 
natures. Based on the previous theories of BI, II, and SI given by the 
formerly referenced studies and on the similarity among BI, II and SI 
discussed above, the basic nature of SI, can be summarized as follows: 

i) For an inviscid fluid, the instability simply depends on the 

2 2 
F* S 

external constraints, i.e., the stability tensor n = ( ) 

S 2 N 2 

which determines the orientation and strength of the 
instability, and has no scale selection. The sufficient 
condition for instability is 

min (N 2 ,F 2 ) < 0, i.e., N 2 F 2 - S 4 < 0 , (2.9) 

2 2 

where the unusual case of both N <0 and F < 0 is not 
included in the second formula. The most unstable slope 
angle is given by (2.4) or (2.7), which with (2.9) leads to 

|k v | < | tan a| < |k Q | . (2.10) 
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Here k v = -F 2 /S 2 is the slope of iso-momentum (i.e., 
fx + V = const) surfaces and also the slope of the 
background absolute vorticity vector, while ke = -S^/N^ i s 
the slope of background isentropic surfaces, 
ii) For a viscous fluid, the instability results from a 
trade-off between two energy transport regimes - the 
conductive regime (no perturbation) and circulation 
regime. When the latter becomes more effective, 
instability occurs and the circulation depends on two 
other competitive factors: the external effect of the 

stability tensor and internal effect of the stress tensor 
due to viscosity and pressure gradient. In this case, a 
scale selection mechanism exists and depends on the 
competition of the above two factors. Since the Prandtl 
number = 1 here, by the already known theory (Chandrasekhar, 
1961 and Veronis, 1966) which asserts that the primary 
bifurcation point must be a steady one unless both rotation 
(measured by the Coriolis parameter) is strong enough and 
the Prandtl number is small enough (much smaller than unity), 
we may expect that the viscous SI motion with certain boundary 
condition will first occur as a steady circulation. 

3. Parcel Dynamics and Energetics 

We shall for simplicity omit the effect of viscosity and thermal 
conductivity. In this case (v = 0), (2.1c-d) give 

4r (v + M) = 0 , (3.1a) 


7 


d_ 

dt 


(3.1b) 


(9 + B J 2 77 < 9 + § ln0 > = 0 » 

o dt o 


where M - fx + V is the absolute momentum of the basic state. If a 
parcel (here a parcel in x = (x,z) plane is a tube along y direction in 
the 3-dimensional space) initially with v = 9 = 0 at moves to 5, 
then the integration of (3.1a-b) following the parcel gives 

v = M(x <> ) - M(x) ,. (3.2a) 



(3.2b) 


(3.3a) 



(3.3b) 


where (’) =— . Since the pressure perturbations only redistribute the 
perturbation energy (cf. Chapter IV, Xu, 1984), we may ignore the effect of 
pressure perturbations and thus follow the methods of Emanuel (1983), to 
obtain physical insights about SI. 
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Alternatively, Eq. (2.5) with v = 0 suggests that the above two step 
analysis, i.e., (3.1) - (3.3), can be simplified by locally introducing a 
"thermal-inertial potential" 


A = f(fx + V) cox a - B 0 sin ct 


for the basic state. It is easy to prove that A = A(x) is locally independent 

A A /S 

of z, i.e., an iso-A line is locally parallel to z coordinate (see Fig. 3.3), 

.A A A A. 

and that 3A/3x = F^. The quantity a defined in (2.3) is the perturbation 
of A., Thus, similarly to (3.2) - (3.3), we have 


a = A(x°) - A(x) , 


and 


'i 



If we consider a, A, x and £ as 6, © 0 , z and S respectively, then we have 
pure BI instability. 

By considering the energy exchange between a displaced parcel and the 

environment, some useful insights into the nature of SI can be gained. For 

2 2 
F S 

a locally uniform basic state, i.e., n = ( ) is constant (3.2a~b) can 

, ... S N 

be rewritten as 


h = -n • l 


(3.4) 
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where h = (fv,9) and £ = (5,0 = |£|(cos a, sin a) is the parcel displacement 

A A 

in the cross section (x,z) or (x,z). Similarly (3.3a~b) can be rewritten as 


d 


2 


dt 


2 


4 - -IT • 


£ - a Vp 
o r 


(3.5) 


Eq. (3.5) indicates that a fluid parcel with displacement i is 
accelerated by two forces: (a) related to the stability tensor of 
the basic flow -H • £, and (b) due to the perturbation pressure 
gradient. Thus, the energy absorbed (or released) by a fluid parcel 
from (or to) the basic flow is 


-A = -P Si 
o 


d* = - 


A f d(F 2 ? 2 + 2S 2 C5 + N 2 C 2 ) 



(JL) 


2 


or 


= - (F^£^ + N 2 £^) in (x,z) coordinates of (2.3), (3.6) 


where ^ = m^(a) is given in (2.8). Eq. (3.6) indicates that (a) -A 
depends on the parcel’s displacement £ only and not on the particular 
trajectory; (b) for fixed |4|, when i is along (or normal to) the 
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most unstable slope of (2*4), i.e*, tan~l(C/£) = a = a (or tan“l(C/5) 

-ala), -A reaches a maximum (or minimum) where ^ = min <*>2 = min 

(F^,n 2) ( or ^2 = max oj 2 = max(F^ , ) ) . 

a 

Furthermore the instantaneous power gained (or lost) by a fluid parcel 
with velocity v at £ from (or to) the basic flow is 

-W - -P y • n •• 1 - - P 0 “» 2 (a,e)|v| 1*1 , |a - 0| < - (3.7) 

where 

aj2(a,P) = cos a cos 8 + sin (a + 0) + sin a sin 3, 

a and 8 are the direction angles of £ and v, respectively* Here the constraint 

| a - 8 1 < -r is required by the fact that for any growing disturbance the 
• d» ; 

angle between v = — and £ , i.e*, |<* - 8|, must acute* In this case, 1 
it is also easy to prove that for any fixed |y| and \i\ -W reaches the 

A A 

maximum (minimum) when a=0=a(a=81a). 

4. Moist SI and CSI 

As mentioned before, the parcel dynamics associated with (3.2a~b) and 
(3.3a~b) does not necessarily require the thermal wind relation between 0 O and 
V. Thus by replacing with 0 W the same analyses and results for the parcel 
dynamics will be exactly applicable to moist SI. However, to estimate the 
most unstable slope angle for moist SI, we need to analyze the moist parcel 
energetics as follows. 

For a moist atmosphere, if the basic moisture field is also only a 
function of (x,z), then in addition to (2.2) there are another two basic 
frequencies N 0 and S Q defined as 



o 


30 



o 


30 

g w 

0 3x 
o 


(4.1) 
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where O w is the wet bulb potential temperature of the basic state. 

If the atmosphere is saturated everywhere, i.e., there is enough con- 
densed liquid (or solid) water in downdrafts, then and in (2. Id) 

2 2 2 

should be replaced by N and S respectively while S in (2.1c) is 

o o 

unchanged. Since does not usually satisfy the thermal wind 

2 2 
K S 

relation, the moist stability tensor ( ~ ~) is generally nonsymmetric 

So No 

and cannot be diagonalized by an orthogonal transformation such as 
(2.3). However, the moist stability tensor always can be transformed 
into either a diagonal or Jordan form by a proper affine transforma- 
tion. The V-derivative terms j especially the nonlinear terms, will 
generally become more complicated. There is no general way to 
transform the moist governing equations into a simple form like (2.5). 
Another difficulty with moist SI is that the integrand £ • H • d£ in 
(3.6) is no longer a total differential if H is nonsymmetric, i.e., 

H * n T . This prevents us from following the parcel energetic analysis 
of (3.6). However, notice that (3.4) and (3.5) are still valid for 
nonsymmetric constant H. Thus we can analyze moist SI with the 
instantaneous parcel energetics of (3.7). 

For nonsymmetric constant II, the o>2(a,g) in (3.7) becomes 

2 2 2 2 
(u (a,B) = F cos oc cos 3 + S cos a sin 3 + S sin a cos 3 

o 

+ sin a sin 3 , |a - &| < y (4.2) 
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It is easy to prove that w (a, 3) reaches its minimum at a = and 

A A A 

6=3, where a and 3 are the roots of 
o o o 


tan (a - 3) 
tan (a + 3) 


2 2 2 2 
(S Z - S Z )/(F Z + N ) , 

o o 


2 2 2 2 
(S* + S )/(F - N Z ) 

o o 


(4.3a) 

(4.3b) 


2 2 * 

For S * S , and thus a * 3 , maximum power (corresponding to 

o o o 

2 

min w (a, 3)) can only be gained by the parcels with velocity 
a,e 

v = (cos sin 0 Q )|y| and displacement 4 = (cos sin <* o )|j|« 

A A 

Since the velocity v generally (except for a = 6 in the case of 

o o 

2 2 

S = S q ) deviates the parcel from ct^ when the maximum power is reached, 

this maximum power can only be gained instantly by the parcel. In 

order to keep the parcel position angle even temporarily constant, the 

parcel velocity must be in the same direction as a, i.e., a = 3. In 

this case the most unstable slope angle a is given by (4.3b) with 

2 2 2 2 

a = 3, or equivalently by (2.4) with + S and instead of 2S and 

2 a 2 - 2 2 2 

N respectively. Thus, N and F are given by (2.6) with + S and 

2 2 2 
N instead of 2S and N respectively, 
o 

We can also show that above analysis and results are consistent 
with following inviscid linear analysis. For a free mode where the 

stream function ~ e ^^ a)t ” mx nz ^ tt is eas y to show that 




cos 2 a + (s 2 + S 2 ) cos ol sine i + N 2 sin 2 c 
o o 


(4.4) 
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where a is the slope angle of stream line, i.e., 


tana * 


♦ = const 


m 

n 


2 2 2 2 2 

Eq. (4.4) is just (2.8) with S + S q and replacing 2S and N respectively. 

Thus, the linear most unstable (or least stable) slope angle and the 

2 

corresponding value of min w are consistent with the above parcel energetics 
result. 

Now we can consider the CSI problem qualitatively (cf.. Fig. 2). In 

dry subsidence regions the downward motion tends to be along the least (dry) 

- 2 2 2 2 
stable slope tana^ which is given by (2.4) with N = and S = where 

the basic state is dry symmetric stable, i.e.. 


2 2 4 ■ 

N 1 F " s i > 0 


(4.5) 


Condition (4.5) implies that the downward motion is passive and must be forced 
in compensation to the moist updrafts. Thus, the basic state in the moist 
updrafts must satisfy moist SI condition, i.e.. 


2 2 2 2 2 

N Z F - (S7 + S ) /4 < 0 
o 1 o 


(4.6) 
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Here the updrafts tend to slant along the most (moist) unstable slope 

- 2 2 2 

tana which is given by (2.4) or (2.7) with S" + S and N replacing 
o 1 o o r 

2 2 

2S and N respect ivley, while the dry subsidence tends to follow the 
slope tana^, thus the averaged slope (angle 4>) is between the two, 
i.e., (see Fig. 2) 

hj < 1^1 < l“ 0 l • (4-7) 

Furthermore, under conditions (4.6) and (4.7), it is easy to prove 


that 


0 < l^j <| \\ <|4> v | < |cT| < { 


4> | < ^ , as N 2 > 0 , 

o' 2 ’ o 


-T<|4>|<Tf,asN 2 <0 , 

2 1 o' o 


(4.8) 


where $]_ and 4> 0 are the slope angles of dry and moist isentropic liens 
(in x plane), <t> v is the slope angle of iso-momentum line. Thus, if <J> equals 
to the slope angle of a moist updraft, then the estimation of (4.7) can 
be improved as 


4> v I < hi < ho 


(4.9) 


because the unstable slope for a rising moist parcel must be steeper than the 
iso-momentum lines (cf. Section 3). Obviously, for a given scale height H, 
(4.7) or (4.9) renders a broad mesoscale range of permissible values for the 
intrinsic horizontal scale L = H tan ^ for CSI circulations (see Fig. 2). 

5. Concluding Remarks 

Via a proper rotational coordinate transformation in a transverse 
section of the basic flow, we have shown that the nonlinear viscous 
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equations with a SI (Symmetric Instability) basic state can be locally 
transformed into equations with a BI (Buoyancy or Convective Instability) or 
II (Inertial Instability) basic state. This indicates the common nature of 
SI, BI and II with respect to their local dynamics. However, it is worth to 
note that for BI the hori zonal stability, due solely to Coriolis force, is 
generally much weaker than the vertical instability. For SI a stable 

a 

stratification can render the ‘’horizontal" (along x) stability much stronger 

A 

than the '’vertical” (along z) instability. This can make for a large 
difference between BI and SI (especially when the Prandtl number + 1, cf. 
Miller, 1984). 

The similarity between SI, BI and II is conducive to a better 
understanding of various SI problems and suggests qualitative results 
for these problems without even solving them analytically. However, 
this analogy does not mean that a SI problem can always be treated as 
easily as a BI or II problem. The additional difficulties with SI 
problems often arise due to the following factors: a) when the basic 
state is not uniform, i.e., the stability tensor is a function of 
position, then the SI problem generally cannot be transformed into a BI 
or II problem; (b) even if the basic state is uniform and thus the SI 
problem can be globally transformed into a BI problem, the boundary 
will usually become more irregular and the problem may still not be 
amenable to analytical treatment. Thus, as far as the global dynamics 
are concerned, SI is not generaly equivalent to BI and II. 

Furthermore, when the instability is conditional, the stratification 
and thus stability tensor will change discontinuously at the interfaces 
between the moist updrafts and dry subsidences. Thus, the differences 
between CSI and CBI (Conditional BI) will be even n*ore than those cited 
above. For CBI, the stability tensor does not have a sudden rotation at 
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the interfaces as it does for CSI. Since the most unstable direction for 


the moist stability tensor is steeper than the least stable direction for 
the dry stability tensor, the CSI circulations will be slantwise between 
these two directions* If the slope of the moist updraft is considered as 
a measure of the slope of a CSI circulation, then it must be also steeper 
than the iso-momentum surfaces of the basic state. Associated with the 
slope of a CSI circulation, the intrinsic horizontal scale may cover a 
broad mesoscale range. These qualitative results are confirmed by linear 
CSI solutions (Chapter III, Xu, 1984). 
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Figure Captions 


Figure 1 Schematic of uniform basic state and rotational coordinate 
transformation R** 

Figure 2 Schematic comparison of the slope angles: 

A 

a - most unstable slope angle in moist basic state, 
a° - least stable slope angle in dry basic state, 

<t> — slope angle of a CSI circulation, 

H - scale height, 

L - intrinsic horizontal scale, 
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*These papers were supported under an associated contract sent directly to 
NASA MSFC which was originally meant to be part of the original contract 
NAS8-33797. 
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Papers submitted or to be submitted are: 

Q. Xu: Conditional Symmetric Instability and Mesoscale Rainbands, 

Quart. J. Roy. Meteor. Soc. » 

Q. Xu: Generalized Energetics for Linear and Nonlinear Symmetric 

Instability, J. Atmos. Sci. . 

C.A. Staver, C.-L. Shie, J.H.E. Clark: Ageostrophic Motions as 

Diagnosed by a Semi-Geos trophic Model and their Influence on 
Cyclone Development, Mon. Wea. Rev. . 

II. Major Findings 

1. We carried out a number of diagnostic studies of the effect of latent 
heat release on a March, 1978 storm. In the M.S. thesis by Staver, we found 
that the large-scale latent heating had a small influence on enhancing the 
growth of the storm. The pattern of quasi-geostrophic vertical motion was 
enhanced in the vicinity of the cold front by the latent heating. We 
estimated the areas of heating from sample measurements and satellite infrared 
and visible imagery. 

2. In the ensuing M.S. thesis by Shie and the paper by Staver, Shie and 
Clark, we extended our diagnostic scheme to a semi-geos trophic formalism. In 
this scheme, advection of heat and vorticity by the ageostrophic as well as 
geostrophic part of the horizontal wind field was accounted for. In the 
thesis, the equations were solved in transformed geostrophic coordinates but 
this proved to be inaccurate because of the many smoothing and interpolation 
operations involved in the coordinate change. In the paper we develop an 
accurate iterative scheme which can be used in physical space. The method 
involves an initial guess of the ageostrophic winds using the quasi- 
geostrophic formalism and then a reevaluation of the advection terms with the 
total wind. Convergence was found to be rapid. We noted only small 
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differences in the large-scale vertical motion field with and without latent 
heating between the semi- and quasi-geos trophic diagnostics. However, with 
the semi-geostrophic total winds we are able to evaluate the tendency of the 
synoptic-scale motion field to produce conditions favorable for symmetric 
instability. This would manifest itself in the form of mesoscale rainbands 
oriented parallel to the low-level vertical wind shear. We noted that within 
a period of less than 12 hours, differential, horizontal and vertical 
advection by the total wind field would render the atmosphere convectively 
unstable along iso-momentum surfaces in the vicinity of the warm front in the 
March, 1978 storm. We thus concluded that l.atent heat release on the 
convective and mesoscale probably plays an important role in the development 
of this winter storm. 

3. We proceeded to investigate some models to account for the 
interaction of convective scale motions and the larger scale motion field in 
which it is embedded. In the paper by Xu and Clark, we considered a wave-CISK 
model with the important feature that the convective motions would not respond 
instantaneously to large scale forcing. Rather they would be triggered by the 
forcing but then would go through a life cycle that depended on 

a. the intrinsic mean cloud lifetime, 

b. the amount of available convective potential energy in the 
environment , and 

c. eddy-dissipation. 

We were able to produce wave-like disturbances which strongly resembled 
traveling waves associated with convection in the vicinity of cold and warm 
fronts and in warm sectors where the background shear is weak. 

We have also carried out a number of theoretical studies of conditional 
symmetric instability (CS1). In the enclosed paper by Xu and Clark, we 
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develop an analytical technique for solving the linear problem with and with- 
out latent heat release. In ensuing papers we investigate some implications 
of nonlinearities and viscosity on the growth of CSI disturbances in an 
unstable environment • 
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